Origin of the multiferroic spiral spin-order in the i^MnOs perovskites 
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The origin of the spiral spin-order in perovskite multiferroic manganites JJMnOs (RE = Tb or 
Dy) is here investigated using a two eg-orbitals double-exchange model. Our main result is that 
the experimentally observed spiral phase can be stabilized by introducing a relatively weak next- 
nearest-neighbor superexchange coupling (~ 10% of the nearest- neighbor superexchange) . Moreover, 
the Jahn- Teller lattice distortion is also shown to be essential to obtain a realistic spiral period. 
Supporting our conclusions, the generic phase diagram of undoped perovskite manganites is obtained 
using Monte Carlo simulations, showing phase transitions from the A-type antiferromagnet, to the 
spiral phase, and finally to the E-type antiferromagnet, with decreasing size of the R ions. These 
results are qualitatively explained by the enhanced relative intensity of the superexchanges. 

PACS numbers: 75.80.-|-q, 75.47.Lx, 75.30.Kz, 64.70.Rh 
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I. INTRODUCTION 

Perovskite manganites, one of the main families of 
strongly correlated electronic materials, have drawn 
much attention since the discovery of the colossal mag- 
netoresistance in the last decade. The strong coupling 
between spin, charge, orbital, and lattice degrees of free- 
dom gives rise to many competing phases in mangan- 
ites, with rich physical properties. Theoretically, the 
physics of manganites appears to be qualitatively under- 
stood within the framework of the double-exchange (DE) 
model, including the superexchange (SE) between the i2g 
spins, and the Jahn- Teller (JT) interactions.— 

Recently, the multiferroic materials, in which the fer- 
roelectric (FE) and magnetic orders coexist and are inti- 
mately coupled, have become a subject of much attention 
due to their technological relevance and fundamental sci- 
ence challenges,^'? The discovery of multiferroicity in un- 
doped manganites with small size R (rare-earth) cations, 
e.g. TbMnOa and DyMnOa, have made the family of 
manganites even more fascinatingfii^ With decreasing 
temperature, these multiferroic materials iJMnOs first 
transform from a paramagnetic (PM) state to a coUinear 
spin sinusoidal incommensurate state at Tn (~ 40 — 50 
K), and then to a spiral spin state with a locked period 
at Tiock (~ 15 - 30 K)^^^ Tjock is also the ferroelec- 
tric (FE) critical temperature (Tc), suggesting a strong 
magnetoelectric interaction. In the coUinear spin sinu- 
soidal incommensurate state (Tlock = T^ < T < Tn) of 
Tbi_a;Dya;Mn03, all spins approximately point into the h 
direction and their amplitudes are modulated by the Mn 
positions in a sinusoidal form. The wavelength of this 
sinusoidal modulation, which is incommensurate to the 
lattice spacing, changes monotonously with decreasing 
temperature until Tiock is reached. Below Tiock, the spins 
form a coplanar (6 — c) spiral order, and its wavelength 
is nearly independent of temperature^ 



Such a spiral spin-order (SSO) driven improper fer- 
roelectricity is also observed experimentally in other 
transition-metal oxides besides undoped manganites, 
such as NiaVsOg,^ CuFeOafi^ CuO^ and others J^ How- 
ever, the theoretical understanding of their properties is 
still in its early stages. For example, the origin of the FE 
polarization remains under much debate. Phenomenolog- 
ically, the SSO breaks the spatial inversion symmetry and 
allows the emergence of a spontaneous polarization.-^ Mi- 
croscopically, two possible mechanisms for the FE polar- 
ization were proposed: (a) a pure electronic contribution 
driven by the spin-orbital couplingfi^ii^ and (b) a cation 
displacement driven by the Dzyaloshinskii-Moriya (DM) 
interactionJ^iil Interestingly, both mechanisms give the 
same behavior: a FE polarization P c>c — e^.j x (Si x Sj), 
with e^j being the unit vector connecting the nearest- 
neighbor (NN) spins Si and Sj . Although recent studies 
suggest the cation displacement to be the dominant one 
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further work is necessary to clarify 



the nature of the FE polarization in these materials. 

In addition, it is important to remark that the origin of 
the SSO is also a puzzle. A direct but phenomenological 
route to generate a SSO phase is via the magnetic frustra- 
tion between NN ferromagnetic (FM) and next-nearest- 
neighbor (NNN) antiferromagnetic (AFM) interactions)^ 
e.g. via a J1-J2- J3 model with classical spins, where Ji 
is the NN superexchange while J2 (J3) is the NNN su- 
perexchange along the b (a) directionsi^ In this model, 
a strong J2 (compared with Ji ) is required to generate a 
realistic short-wavelength spiral order. In fact, the min- 
imum J2/IJ1I for the spiral phase is 0.5 and to reach a 
realistic wavevector q, then J2/I Ji| should be larger, e.g. 
0.8 ~ 1.^^ However, in first-principles calculations^^ the 
value of J2 was found to be only 0.56|Ji|, which is not 
sufficient to induce the experimentally observed short- 
length (L = 6) spiral order. Also, the coupling J3 was 
found to be 0.375| Ji| and in addition it is AFM, opposite 



to the weak FM coupling (O.OlJi) used in the previous 
model studyf^i Furthermore, while providing interesting 
information for the spiral regime, it should be noted that 
the classical spin model is not suitable to explore in a 
single framework the many phases (A-AFM, spiral, E- 
AFM) found in RMnOa manganites, since it can not dis- 
tinguish the E-AFM order from the q = 1/4 spiral order 
(see Fig. 1) since they have the same energy. In sum- 
mary, while the classical spin model with NNN magnetic 
frustration provides a natural starting point to describe 
phenomenologically the SSO phase, this simple model is 
not sufficient to deeply understand the microscopic origin 
of the SSO in perovskite manganites. 

An alternative route to obtain the SSO phase is to in- 
corporate the DM interaction (ex Si x Sj ) into the double- 
exchange framework. ^^'^'^ However, the required intensity 
of the DM interaction is two orders of magnitude higher 
than expected.^'' Thus, it is fair to express that the real 
driving force for SSO in manganites remains a mystery. 
Note that a more fundamental theory to explain the ori- 
gin of the SSO phase should also explain the general 
phase diagram of i?Mn03: with decreasing R size, the 
ground state of i?Mn03 changes from the A-type AFM 
(A) state, as in the case of LaMnOa, to the spiral (S) 
state, as it occurs in TbMnOs, and finally to the E-type 
AFM (E) phase, which is stabilized in HoMnOs. This 
property is referred below as the "A-S-E transition" i^^ 
For all these phases, the order of the Mn'^+ spins can 
be characterized by a particularly propagation vector (0, 
^Mn, 1) (in the orthorhombic Pbnm cell notation, see 
Fig. 1), indicating an AFM coupling along the c axis, 
and FM coupling along the a axis. ^Mn is for the A 
phase and 0.5 for the E phase. However, for convenience 
when using a square lattice as in our simulations, q in 
the following will be defined along the x/y directions (see 
Fig. 1) and it equals half of the gMn used in experimental 
papers. 
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FIG. 1: (Color online) (a) Sketch of the crystal structure 
(a — b plane) of RMnOi. Two types of coordinate axes (a — 
b and x — y) are shown, (b) Illustration of the two kinds 
of distortions discussed in the text: GdFeOa-type (oxygen 
moves perpendicular to the Mn-O-Mn bond) and Jahn- Teller 
type (oxygen moves along Mn-O-Mn bond), (c) Two Jahn- 
Teller distortion modes: Q2 and Q3. (d) Sketch of the spin 
patterns (in the a—b plane) corresponding to seven (5*4 and 5*6 
belong to the same group) candidate phases considered here 
in the zero-temperature variational method. A, G, CE, C, 
and E are standard notations for well-known AFM phases in 
manganites. Here, the ferromagnetic chains in the C/CE/E- 
AFM phases are within the a — b plane. D denotes the TTii 
dimer phase, that is considered in our studies although it has 
not been observed in real manganiteS;-- In the spiral phases 
SL, the length period is denoted by L (along x/y). The L = 4 
case has the same energy as the E-AFM phase in the classical 
spin model context. 



II. MODEL AND METHODS 

In this study, we will analyze the properties of the two 
Cg-orbitals DE model for manganites defined on a two- 
dimensional (a-b plane) lattice to try to unveil the ori- 
gin of the SSO phase in i?Mn03. We will start from a 
pure DE model, and then incorporate other interactions 
one-by-one to clarify their respective roles. The primary 
Hamiltonian, considering only the DE and NN SE inter- 
actions, reads as: 
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where the first term is the standard DE kinetic en- 
ergy of the Cg electrons. The DE hopping amplitudes 
f^^ are orbital- and direction-dependent. In particular, 
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to where the superscript 1 (2) denotes the eg 



orbital dx^-y2 (d^z'^-r^), and to (0.2 ~ 0.3 eV) is taken 
as the energy uniti The infinite Hund coupling used 
here generates a Berry phase flij — cos(6',j/2) cos{9j/2) + 
sm(0i/2) sin(9j/2) exp[—i(ipi — (pj)], where 9 and (p are 
the angles defining the t2g spins S in spherical coordi- 
nates. The Cg kinetic (DE) energy is obtained by exactly 
diagonalizing the fermionic sector, using library subrou- 
tines, once a t2g classical spins configuration is generated. 
The second term in the Hamiltonian is the usual AFM 
SE couphng between NN t2g spins. 

The numerical methods used in this manuscript are of 
two varieties: the variational method employed at zero 
temperature and the Monte Carlo (MC) simulation at 
finite temperature. For the case of the zero tempera- 
ture variational method, the total energies (per site) of 
several candidate phases (for an infinite size lattice) are 
calculated to determine which is the most likely ground 
state. The candidate phases are sketched in Fig. 1, and 
almost all of the typical spin order patterns discussed 
in manganites have been included. In addition, spirals 



with wavevectors q (from to 1/4) are also taken into 
account. For the MC simulation, the classical i2g spins 
evolve in MC time following the standard Metropolis al- 
gorithm. The MC simulation is independent of the zero- 
temperature variational method, namely the spin order 
obtained via MC simulations can be one of the candi- 
date phases, but the outcome is unbiased and it is not 
limited to be one of those phases. Due to the restriction 
of working on a finite lattice size, a large enough cluster 
should be used in the MC simulation to allow for other 
possible spin orders. Since the MC method used here 
has been discussed extensively in the manganite litera- 
ture, for more details regarding the Hamiltonian and the 
MC technique the readers are referred to Ref. HI 



to 6 > a. For a crude estimation, if all Mn'^+ and 0^~ 
in Fig. 1 are simplified to be coplanar, it is easy to cal- 
culate the bond angles 4>h and 0q from the experimental 
data available for a, 6, and </>. For the case of TbMnOs, 
it was found that h/a « 1.1 and (f> « 145.3° suggest- 
ing (t)b w 115.2° and 0a w 99.5°,^i^ The larger bond 
angle leads to a stronger NNN spin coupling through 
the Mn(l)-0-0-Mn(4) exchange path. Considering the 
angle-dependence of the super-exchange^^ the NNN SE 
J2b coupling between Mn(l) and Mn(4) can be between 
1.8 to 2.2 times the value of the J2a coupling between 
Mn(2) and Mn(3). In addition, the Cg orbital-order can 
also induce this anisotropy^ Thus, the new term in the 
Hamiltonian incorporating the NNN SE coupling is 



III. ZERO TEMPERATURE RESULTS 

For the model studied here, the Cg electronic den- 
sity is always 1 when the chemical potential equals 0. 
Thus, in the "clean limit" two-orbitals DE model with- 
out quenched disorder, the undoped RMjiO^ case is free 
from electronic phase separation tendencies, which usu- 
ally occurs in doped casesi Therefore, it is reasonable to 
consider only homogeneous phases as candidates for the 
ground states. A very important result in this context 
is that using Eq. 1 we have observed that the total en- 
ergy of the spiral phase is always higher than the energy 
of both the A and the E phases, regardless of the value 
of Jaf- In other words, the spiral phase can not be the 
ground state of Eq. 1, in agreement with previous studies 
of undoped manganites.^^ This model has to be modified 
for the spiral state to become the ground state in some 
range of parameters. 



A. The role of the NNN superexchange 

In the rest of the paper, the spin frustrating effects 
of the NNN J2 coupling will be taken into consideration 
to solve this paradox. Here, J2 arises from the super- 
exchange between NNN 3d spins. It should be noted 
that our model is conceptually different from those based 
on classical spins presented in previous publicationsi^i^ 
First, our NN SE is antiferromagnetic while it was fer- 
romagnetic in those previous efforts. Ferromagnetic ten- 
dencies do exist in our model but they originate in the 
DE mechanism of itinerant eg electrons. Second, as 
shown in Fig. 1, considering the substantial distance be- 
tween the NNN spins via the zigzag exchange path Mn- 
0-0-Mn, then J2 should be weaker than Jaf, in con- 
trast to the robust J2 used before^^i^ As indicated in 
Fig. 1, the GdFeOs-type distortion of the oxygen octa- 
hedra, which reduces the NN Mn-O-Mn angle 0, is cru- 
cial for J2, because it shortens the distance between two 
oxygens and increases the Mn-0-O-Mn angles. There- 
fore, J2 is enhanced by the GdFeOa-type distortion,— 
Another important property of J2 is its anisotropy due 
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where 7 runs over a and 6, and [ij] denotes NNN sites. 
Considering together Eqs. 1 and 2, a zero-temperature 
phase diagram can be obtained by comparing the total 
energies of the candidate phases discussed in Fig. 1(d). 
As shown in Fig. 2(a), it is remarkable that now the spiral 
phase becomes the ground state for values of J^h as small 
as 0.017. 

Considering the GdFeOs-type distortion, we calculated 
the 0-dependent Jaf (using the relation Jaf csc sin (0/2) 
described in Ref. |25|) and also the DE kinetic energy (us- 
ing t^'^ ex cos^ 0, i^'^ and i^'^ ex cos^ 0, and f^'^ ex cos0, 
as in Ref. uw of the A and E phases. The latter are 
denoted by K^ and K^, respectively, and are shown in 
Fig. 2(b). Both the DE kinetic energy and Jaf decrease 
with decreasing 0. The curves for Kp^ and K-^ are al- 
most identical and lower than the Jaf curve. Therefore, 
the GdFeOs-type distortion increases Jaf relative to the 
DE kinetic energy. Noting that J2 is enhanced by the 
GdFeOs-type distortion, then the total effect of this dis- 
tortion is the enhancement of both the NN and NNN SE 
couplings, particularly the latter. Therefore, the A-S-E 
phase transition with decreasing R size can now be easily 
understood, and the arrow in Fig. 2(a) indicates a pos- 
sible path that is qualitatively in agreement with experi- 
ments. Below, we will still use the ideal case = 180°, as 
in most previous studies, since the GdFeOs-type distor- 
tion effect can be accounted for by increasing the relative 
values of Jaf/^o and J2/to (see Fig. 2b). 

As shown in Fig. 2(c), at the A to S transition the 
wavevector q changes continuously from to a finite 
value. However, between the spiral and E phases, not 
only q now changes discontinuously, but also their spin 
patterns are incompatible, namely the E-phasc is not 
a special case of the spiral phase. These two different 
phase boundaries may cause two distinct behaviors when 
R cations substitutions occur: (a) if Tb is substituted by 
an A-phase R (e.g. Gd), q will be reduced continuously; 
(b) if Dy is substituted by an E-phase R (e.g. Ho), phase 
separation at some concentrations may be induced. 
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FIG. 2; (Color online) (a) Zero-temperature phase diagram 
of the two-orbitals DE model for iJMnOs. The possible path 
for the A-S-E phase transition is indicated by the arrow. The 
phase diagram is independent of J2a, as long as J 2a < J2h- 
(b) Jaf, Ka and Ke, as a function of the Mn-O-Mn angle, 
normalized to their values ai (j> ~ 180°. (c) Wavevector num- 
ber q for (a). Here, q is defined along the x/y direction and 
equals half of the qMn used in experimental papers. 
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FIG. 3: (Color online) (a) Zero-temperature phase diagram 
of the two-orbitals DE model with JT distortions (A|Q2|=1.5) 
for i^MnOs. The possible values of JAF-J2h for realistic gs in 
TbMnOa and DyMnOa are shown with dashed lines. The 
four asterisks represent the set of couplings (1st and 2nd for 
A, 3rd for S, and 4th for E) that are studied via the MC 
techniques, (b) Same as (a) except for a finite J2a=0.5J2b. 
The four asterisks in (a) remain within the A-S-E region here. 
(c) Wave vector number q for (a). 



B. The role of the Jahn- Teller couplings 

In spite of the success in describing the SSO phase in 
the previous section, we have observed that in Fig. 2(c) 
the wavevector q is still smaller than the experimental 
value. For example, to obtain the q for TbMnOs (0.14 
along the x/y direction), the minimum J2b is 0.073 which 
corresponds to Jaf=0.074. This J^h/ Jaf ratio is too 
large for real manganites, as mentioned before in this 
text. Moreover, the largest q in Fig. 2 can not reach the 
experimental value for DyMnOa (g=0.19). Even if the 
GdFeOs-type distortion is taken into account (0wl44.7° 
for DyMnOa),^ q can not be increased further. 

Interestingly, we have observed that realistic values for 
q can be stabilized if the Jahn- Teller distortion (Fig. 1(b)) 
is incorporated. This JT distortion can be characterized 
by two modes: Q2 and Qs (Fig. 1(c)), which couple with 
the e„ electrons via: 



-ff,TT = ^2^[Q2,tTx.i + Qz,iTz 



(3) 



where A is the spin-phonon coupling coefficient and r is 
the orbital pseudospin operator, given by t^ = c\ch+clca 



and Tz = c\ca — c\^Cht^ For all i^MnOs at low tempera- 
tures, \Q2\ and Qs are uniform with \Q2\ ~ — vSQs- The 
sign of Q2 is staggered, which gives rise to the well-known 
staggered d^y^'2_^2 and d3y2_i.2 orbital order JS Therefore, 
Eq. 3 can be scaled using AIQ2I as the only parameter. 

Considering together Eqs. 1, 2 and 3, we recalculated 
the zero-temperature phase diagram. A typical result 
with A I (32 1 = 1-5 is shown in Fig. 3. Compared with 
Fig. 2(a), now the spiral phase region in Fig. 3(a) is en- 
larged. This spiral phase can even be stable without J^hi 
although in a very narrow Jaf region. In Fig. 3(c), the 
wavevector q of the spiral phase is increased as a whole, 
compared with those in Fig. 2(c). Therefore, the q for 
TbMnOa and DyMnOs can be obtained with a small 
J2b, e.g. the minimum J2b is now 0.002 (Jaf=0.091) 
for TbMnOg, and 0.008 (Jaf=0.09) for DyMnOg, re- 
spectively. In other words, with the JT distortion incor- 
porated, a rather weak NNN SE (< 10% of NN SE) is 
enough to generate the realistic SSO in iJMnOs. In con- 
trast to the case without JT distortions, J2a also plays 
a prominent role here, although it can not change the 
phase boundary between the A-S-E phases. Moreover, 
Fig. 3(b) shows that a robust J2a even induces C and D 
phases in the large J2b region, implying that both J2b 



and J2a in real manganites should be weak since none of 
these two phases have been observed experimentally in 
bulk i?Mn03. 

Besides providing realistic values for q, the JT distor- 
tion can contribute to the insulating nature of -RMnOa, 
which is crucial for the FE polarization. For instance, 
with AIQ2I = 1-5, the calculated gap (not shown) is 
2.59 {^ 0.52 - 0.78 eV) for the spiral phase at the q 
of TbMnOs. The agreement with experimental results 
for the gap (0.489 eV) suggests that the chosen value for 
AIQ2I is reasonable. ■^^ 

To understand the enhancement of q by the JT effect, 
let us focus on the A to S transition. We consider ipi = 
and 9i = 27rq-ri for the <2g spins. After diagonalizing the 
fermionic operators in Eqs. 1-3, the ground state energy 
per site is E ^ K + 2 Jaf cos(A6') + J2h cos(2A6') + J2a, 
where A9 — 2Tiq is the angle between NN spins, and 
the fermionic sector energy is K Ki Kxcos{^). Ex- 
panding E around A0 = 0, we obtain E Ki Eq — (-^ -|- 
JAF + 2J2b)Ae' + (il + ^ + ^)Ae^ + O(A06)^^l^g^g 

Eq is a constant and A6' (hence q) becomes an order 
parameter. The A-S phase boundary is determined by 
—^ + Jaf + 2 J2b = 0, which agrees with the numerical 
results shown in Figs. 2(a) and 3(a-b). From these con- 
siderations, it is clear that the competition between Kp^^ 
Jaf, and J2b is crucial to obtain the spiral phase using 
realistic interactions, distinguishing our model from the 
previously used effective Ji — J2 spin models. Further- 
more, the reduction of —Ka by the JT distortion lowers 
the threshold of (Jaf, J2b) for a spiral phase which, as a 
consequence, leads to an enhancement of q. 



IV. FINITE TEMPERATURE RESULTS 

The above described phase diagrams were obtained at 
zero-temperature comparing the energies of several can- 
didate phases. In principle, this procedure can not rule 
out the possibility of other unknown phases becoming 
stable. Therefore, it is necessary to check the results 
using MC techniques. For simplicity, the spins are re- 
stricted to be in an easy-plane for the magnetic moments, 
namely we use the X — Y model instead of the Heisen- 
berg model for the i2g spins. This simplification, which 
avoids the introduction of extra adjustable parameters to 
generate a magnetic anisotropy, is reasonable from the 
point of view of the physics of these materials since all 
spins are coplanar for all known low-temperature phases 
in .RMnOa. The JT distortions are fixed as in the pre- 
viously described zero-temperature case. P is obtained 
using the model- free equation — e^j x (S^ x Sj). We use 
a 2-dimensional (2D) L x L {L = 12) lattice with peri- 
odic boundary conditions and this computational study 
is carried out only at four sets of couplings (J2b, Jaf), 
which are denoted by asterisks in Fig. 3(a). The very 
time consuming MC runs were typically carried out with 
6000 steps for thermalization and 4000 steps for mea- 
surements. Considering that all the physics of the four 



coupling sets (asterisks) is not much affected by a weak 
J2a (see Figs. 3(a) and (b)), here J2a is set to be zero for 
simplicity. The q for the 3rd asterisk should be the com- 
mensurate value 1/6 (suitable for Tbo.4iDyo.59Mn03) to 
match the lattice size.- 

For all four pairs of couplings, the presence of sharp 
characteristic peaks in the spin structure factor of our 
MC results (Fig. 4(a)) confirms the stability of the can- 
didate phases at low temperature. In addition. Fig. 4(b) 
shows that the T^ for the FE transition in the spiral phase 
is found to be 0.0065 (~ 16 — 24 K), agreeing with the 
experimental value 22 Ki^ It should be noted that al- 
though the Mermin- Wagner theorem forbids a finite Tc 
(for an infinite lattice size) in the 2D X — Y model, ex- 
tensive MC studies on double-exchange models in the 
past decade have suggested that the characteristic T* in 
2D small size lattices where the order develops (namely 
when the correlation length becomes comparable to the 
lattice size) can still be comparable with the experimen- 
tal critical temperatures, at least qualitatively. In a 2D 
regime of very large correlation lengths, a weak coupling 
in the third direction rapidly establishes long-range order 
with T* as the critical temperature. The real-space cor- 
relation between spins at the maximum distance in our 
lattice C(-j, "I) can also be used to provide estimations 
of the critical transition temperature. 
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FIG. 4: (Color online) (a) Spin structure factor of the four 
MC simulations at T = 0.004, corresponding to the four sets 
of couplings described in the text. AH characteristic peaks are 
prominent, suggesting stable spin orders, (b) Temperature- 
dependent polarization P and C(-|-, ■§■) corresponding to the 
spiral phase with q — 1/6. All values are normalized to their 
saturation values, (c) Sketch of the finite-temperature phase 
diagram. The Tn's for the A and E phases are determined 
from their C(-§, ■|-)-temperature curves. The lines connecting 
Tjv (Tc) are just to guide the eyes. 



Based on these considerations, here a finite- 
temperature phase diagram is sketched (Fig. 4(c)), 
guided by the four couphng sets of Fig. 3(a). It should 
be noted that in principle the coUinear spin sinusoidal 
incommensurate state at Tc < T < Tn is not possible in 
our study, because in our model the spin amplitude 1 5*1 is 
a constant and the period of the phases should be com- 
mensurate with L. Despite this finite-size effect problem, 
our phase diagram nevertheless resembles at least quali- 
tatively the real one, and moreover our rough estimations 
for the Tn of the A and E phases also agree well with ex- 
perimental results 1^ 



including the GdFeOs-type distortion in our study. Note 
that this weak NNN SE interaction (< 10% NN superex- 
change) , while shown here to be crucial in the context of 
the manganite multiferroics, it does not alter the previous 
large body of investigations and conclusions reached via 
similar MC simulations for undoped and doped LaMnOa, 
since in that case the extra NNN SE couplings can be ne- 
glected. 
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